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A theoretical treatment is presented of kinetic equilibria in accretion discs around compact ob- 
jects, for cases where the plasma can be considered as collisionless. The plasma is assumed to be 
axi-symmetric and to be acted on by gravitational and electromagnetic fields; in this paper, the 
particular case is considered where the magnetic field admits a family of toroidal magnetic sur- 
faces, which are locally mutually-nested and closed. It is pointed out that there exist asymptotic 
kinetic equilibria represented by generalized bi-Maxwellian distribution functions and characterized 
by primarily toroidal differential rotation and temperature anisotropy. It is conjectured that kinetic 
equilibria of this type can exist which are able to sustain both toroidal and poloidal electric current 
densities, the latter being produced via finite Larmor-radius effects associated with the tempera- 
ture anisotropy. This leads to the possibility of existence of a new kinetic effect - referred to here 
as a "kinetic dynamo effect" - resulting in the self-generation of toroidal magnetic field even by a 
stationary plasma, without any net radial accretion flow being required. The conditions for these 
equilibria to occur, their basic theoretical features and their physical properties are all discussed in 
detail. 

PACS numbers: 95.30.Qd, 52.30.Cv, 52.25.Xz, 52.55.Dy, 52.25.Dg 



I. INTRODUCTION 



This paper is concerned with dynamical processes in 
astrophysical accretion disc (AD) plasmas and their re- 
lationship with the accretion process. The aim of the 
research programme of which the first part is reported in 
the present paper is to provide a consistent theoretical 
formulation of kinetic theory for AD plasmas, which can 
then be used for investigating their equilibrium proper- 
ties and dynamical evolution. Note that what is meant 
here by the term "equilibrium" is in general a stationary- 
flow solution, which can also include a radial accretion 
velocity. Apart from the intrinsic interest of this study 
for the equilibrium properties of accretion discs, the con- 
clusions reached may have important consequences for 
other applications and for stability analyses of the discs. 
In this paper we consider the particular case where the 
AD plasma contains domains of locally-closed magnetic 
surfaces where there is in fact no local net accretion. We 
then focus on these domains. We do this, both because 
this represents a situation which is of considerable poten- 
tial interest, but also because it leads to some significant 
simplifications of the discussion. In subsequent papers, 
we will proceed to consider more general cases. 
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A. Astrophysical background 



Accretion discs are observed in a wide range of astro- 
physical contexts, from the small-scale regions around 
proto-stars or stars in binary systems to the much larger 
scales associated with the cores of galaxies and Active 
Galactic Nuclei (AGN). Observations tell us that these 
systems contain matter accreting onto a central object, 
losing angular momentum and releasing gravitational 
binding energy. This can give rise to an extremely power- 
ful source of energy generation, causing the matter to be 
in the plasma state and allowing the discs to be detected 
through their radiation emission [H-0]. A particularly 
interesting class of accretion discs consists of those oc- 
curring around black holes in binary systems, which give 
rise to compact X-ray sources. For these, one has both 
a strong gravitational field and also presence of signifi- 
cant magnetic fields which are mainly self-generated by 
the plasma current densities. Despite the information 
available about these systems, mainly provided by obser- 
vations collected over the past forty years and concern- 
ing their macroscopic physical and geometrical properties 
(structure, emission spectrum, etc.), no complete theo- 
retical description of the physical processes involved in 
the generation and evolution of the magnetic fields is yet 
available. While it is widely thought that the magneto- 
rotational instability (MRI) 3] plays a leading role in 
generating an effective viscosity in these discs, more re- 
mains to be done in order to obtain a full understand- 
ing of the dynamics of disc plasmas and the relation of 
this with the accretion process. This requires identifying 
the microphysical phenomena involved in the generation 
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of instabilities and/or turbulence which may represent 
a plausible source for the effective viscosity which, in 
turn, is then related to the accretion rates [H, 0] ■ There 
is a lot of observational evidence which cannot yet be 
explained or fully understood within the framework of 
existing theoretical descriptions and many fundamental 
questions still remain to be answered [3, Q . 



ature anisotropics, and is essential for making a complete 
study of the kinetic instabilities which can play a key 
role in causing the accretion process [l6l - [T8j . We note 
here that, although there are in principle several possi- 
ble physical processes which may explain the appearance 
of temperature anisotropies (see for example |17l . [HI), 
the main reason for their maintenance in a collisionless 
and non-turbulent plasma may simply be the lack of any 
efficient mechanism for temperature isotropization. 



B. Motivations for a kinetic theory 

Historically, most theoretical and numerical investi- 
gations of accretion discs have been made in the con- 
text of hydrodynamics (HD) or magneto-hydro-dynamics 
(MHD) (fluid approaches) dHJ, lO-flllj. Treating the 
medium as a fluid allows one to capture the basic large- 
scale properties of the disc structure and evolution. An 
interesting development within this context has been the 
work of Coppi [12] and Coppi and Rousseau [l]| who 
showed that stationary magnetic configurations in AD 
plasmas, for both low and high magnetic energy densities, 
can exhibit complex magnetic structures characterized 
locally by plasma rings with closed nested magnetic sur- 
faces. However, even the most sophisticated fluid models 
are still not able to give a good explanation for all of 
the complexity of the phenomena arising in these sys- 
tems. While fluid descriptions are often useful, it is well 
known that, at a fundamental level, a correct descrip- 
tion of microscopic and macroscopic plasma dynamics 
should be formulated on the basis of kinetic theory (ki- 
netic approach) [3 [Uj]> and f° r that there is still, re- 
markably, no satisfactory theoretical formulation. Going 
to a kinetic approach can overcome the problem charac- 
teristic of fluid theories of uniquely defining consistent 
closure conditions [ij], [H| > and a kinetic formulation is 
necessarily required, instead of MHD, for correctly de- 
scribing regimes in which the plasma is either collision- 
less or weakly collisional [16H18I ]. In these situations, the 
distribution function describing the AD plasma will be 
different from a Maxwellian, which is instead character- 
istic of highly collisional plasmas for which fluid theories 
properly apply. An interesting example of collisionless 
plasmas, arising in the context of astrophysical accre- 
tion discs around black holes, is that of radiatively in- 
efficient accretion flows pjj ■ Theoretical investiga- 
tions of such systems have suggested that the accreting 
matter consists of a two-temperature plasma, with the 
proton temperature being much higher than the elec- 
tron temperature [H, This in turn implies that the 
timescale for energy exchange by Coulomb collisions be- 
tween electrons and ions must be longer than the other 
characteristic timescales of the system, in particular the 
inflow time. In this correct physical description 

of the phenomenology governing these objects can only 
be provided by a kinetic formulation. Finally, the kinetic 
formalism is more convenient for the inclusion of some 
particular physical effects, including ones due to temper- 



ed. Previous work 

Only a few studies have so far addressed the problem 
of deriving a kinetic formulation of steady-state solutions 
for AD plasmas. 

based on 
2j for lab- 
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theoretical results obtained by Mahajan [2 
oratory plasmas, is a first example going in this direc- 
tion. These authors assumed an equilibrium distribution 
function expressed as an infinite power series in the ratio 
of the drift velocity to the thermal speed (considered as 
the small expansion parameter) such that the zero-order 
term coincides with a homogeneous Maxwellian distri- 
bution. Assuming a prescribed profile for the external 
magnetic field and ignoring the self-generated field, they 
looked for analytic solutions of the Vlasov-Maxwell sys- 
tem for the coefficients of the series (typically truncated 
after the first few orders). However, the assumptions 
made strongly limited the applicability of their model. 

Another approach proposed recently by Cremaschini 
et al. 15], gives an exact solution for the equilibrium 
Kinetic Distribution Function (KDF) of strongly magne- 
tized non-relativistic collisionless plasmas with isotropic 
temperature and purely toroidal flow velocity. The strat- 
egy adopted was similar to that developed by Catto et al. 
|24| for toroidal plasmas, suitably adapted to the context 
of accretion discs. The stationary KDF was expressed in 
terms of the first integrals of motion of the system show- 
ing, for example, that the standard Maxwellian KDF is 
an asymptotic stationary solution only in the limit of a 
strongly magnetized plasma, and that the spatial profiles 
of the fluid fields are fixed by specific kinetic constraints 

In recent years, the kinetic formalism has also been 
used for investigating stability of AD plasmas, particu- 
larly in the collisionless regime and focused on studying 
the role and importance of MRI [181125^27} . The main 
goal of these studies ^6H18l . [26j was to provide and test 
suitable kinetic closure conditions for asymptotically- 
reduced fluid equations (referred to as "kinetic MHD"), 
so as to allow the fluid stability analysis to include 
some of the relevant kinetic effects for collisionless plas- 
mas [HI, [HI, [26[. However, none of them systematically 
treated the issue of kinetic equilibrium, and the under- 
lying unperturbed plasma was usually taken to be de- 
scribed by either a Maxwellian or a bi-Maxwellian KDF. 
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Finally, increasing attention has been paid to the role 
of temperature anisotropy and the related kinetic insta- 
bilities. Some recent numerical studies [l^,|22j have tried 
to include the effects of temperature anisotropy but, al- 
though it is clear that this can give rise to an entirely new 
class of phenomena, all of these estimates rely on fluid 
models in which kinetic effects are included in only an ap- 
proximate way. A kinetic approach is needed rather than 
a fluid one, in order to give a clear and self-consistent pic- 
ture, and this needs to be based on equilibrium solutions 
suitable for accretion discs. 



D. Open problems 

Many problems remain to be addressed and solved re- 
garding the kinetic formulation of AD plasma dynamics. 
Among them, we focus on the following: 

f . The construction of a kinetic theory for AD plas- 
mas within the framework of the Vlasov-Maxwell 
description, and the investigation of their kinetic 
equilibrium properties. 

2. The inclusion of finite Larmor-radius (FLR) effects 
in the MHD equations. For magnetized plasmas, 
this can be achieved by making a kinetic treatment 
and representing the KDF in terms of gyrokinetic 
variables (Bernstein and Catto (281431} ). The gy- 
rokinetic formalism provides a simplified descrip- 
tion of the dynamics of charged particles in the 
presence of magnetic fields, thanks to the symme- 
try of the Larmor gyratory motion of the particles 
around the magnetic field lines. Therefore kinetic 
and gyrokinetic theory are both fundamental tools 
for treating FLR effects in a consistent way. 

3. The determination of suitable kinetic closure con- 
ditions to be used in the fluid description of the 
discs. These should include the kinetic effects of 
the plasma dynamics in a consistent way. 

4. Extension of the known solutions to more general 
contexts, with the inclusion of important effects 
such as temperature anisotropy. 

5. Development of a kinetic theory for stability anal- 
ysis of AD plasmas. As already mentioned, this 
could throw further light on the physical mecha- 
nism giving rise to the effective viscosity and the 
related accretion processes. This is particularly in- 
teresting for collisionless plasmas with temperature 
anisotropy, since only kinetic theory could be able 
to explain how instabilities can originate and grow 
to restore the isotropic properties of the plasma. 



E. Goals of the paper 

The aim of this paper is to extend the investigation 
of kinetic equilibria developed in an earlier paper [l5j 
to a more general class of solutions. In particular, we 
pose here the problem of constructing analytic solutions 
for exact kinetic and gyrokinetic axi- symmetric gravita- 
tional equilibria (see definition below) in accretion discs 
around compact objects. The solution presented is ap- 
plicable to collisionless magnetized plasmas with temper- 
ature anisotropy and mainly toroidal flow velocity. The 
kinetic treatment of the gravitational equilibria neces- 
sarily requires that the KDF is itself a stationary solu- 
tion of the relevant kinetic equations. Ignoring possible 
weakly-dissipative effects, we shall assume - in particu- 
lar - that the KDF and the electromagnetic (EM) fields 
associated with the plasma obey the system of Vlasov- 
Maxwell equations. The only restriction on the form of 
the KDF, besides assuming its strict positivity and it be- 
ing suitably smooth in the relevant phase-space, is due 
to the requirement that it must be a function only of the 
independent first integrals of the motion or the adiabatic 
invariants for the system. 

The paper is organized as follows. In Section 2 we 
discuss the conditions for the existence of a kinetic equi- 
librium and we introduce the basic assumptions for the 
formulation of the kinetic theory. Section 3 deals with 
the first integrals of motion and the gyrokinetic adia- 
batic invariants of the system. In Section 4 we construct 
the equilibrium KDF for AD plasmas with non-isotropic 
temperature, giving also a useful asymptotic expansion 
for this in the limit of strong magnetic fields. Section 5 
deals with calculation of the fluid moments of the sta- 
tionary KDF. The limit of isotropic temperature is then 
investigated in Section 6, while Section 7 is devoted to 
analyzing the Maxwell equations. Finally, Section 8 con- 
tains conclusions and a summary of the main results. 



II. KINETIC THEORY FOR ACCRETION DISC 
PLASMAS: BASIC ASSUMPTIONS 

We first discuss what is meant by asymptotic kinetic 
equilibria in the present study and what are the physical 
conditions under which they can be realized. 

An asymptotic kinetic equilibrium must be one ob- 
tained within the context of kinetic theory and must be 
described by the stationary Vlasov-Maxwell equations. 
This means that the generic plasma KDF, / s , must be a 
solution of the stationary Vlasov equation, as will be the 
case if f s is expressed in terms of exact first integrals of 
the motion or adiabatic invariants of the system, which 
in turn implies that f s for each species must be an ex- 
act first integral of the motion or an adiabatic invariant. 
The stationarity condition means that the equilibrium 
KDF cannot depend explicitly on time, although in prin- 
ciple it could contain an implicit time dependence via 
its fluid moments (in which case the kinetic equilibrium 
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does not correspond to a fluid equilibrium and there are 
non-stationary fluid fields). 

In the following we shall take the AD plasma to be: a) 
non-relativistic, in the sense that it has non-relativistic 
species flow velocities, that the gravitational field can 
be treated within the classical Newtonian theory, and 
that the non-relativistic Vlasov kinetic equation is used 
as the dynamical equation for the KDF; b) collisionless, 
so that the mean free path of the plasma particles is much 
longer than the largest characteristic scale length of the 
plasma; c) axi- symmetric, so that the relevant dynamical 
variables characterizing the plasma (e.g., the fluid fields) 
are independent of the toroidal angle (p, when referred to 
a set of cylindrical coordinates (R,ip,z); d) acted on by 
both gravitational and EM fields. 

Also, we will focus on the situation where the equilib- 
rium magnetic field B admits, at least locally, a family 
of nested axi-symmetric closed toroidal magnetic surfaces 
{ip{r)} = {V-'(r) = const.}, where ifj denotes the poloidal 
magnetic flux of B (see [H, EH for a proof of the pos- 
sible existence of such configurations in the context of 
astrophysical accretion discs; see also [2 EH for further 
discussions in this regard and Fig.l for a schematic view 
of such a configuration). In this situation, a set of mag- 
netic coordinates (tp, ip, ■&) can be defined locally, where d 
is a curvilinear angle-like coordinate on the magnetic sur- 
faces ip( r ) = const. Each relevant physical quantity A(r) 
can then be expressed as a function of these magnetic 
coordinates, i.e. A(r) — A(ip, •&) , where the <p depen- 
dence has been suppressed due to the axi-symmetry. It 
follows that it is always possible to write the following de- 
composition: A = A^ + {A} , where the oscillatory part 
A~ = A — (A) contains the ^-dependencies and (A) is the 
if)— surface average of the function A(r) defined on a flux 
surface ip(r) = const, as {A} = £ _1 § d&A{r)/ |B • W| , 
with £ denoting £ = §d&/ |B ■ W|. 

For definiteness, we shall consider here a plasma con- 
sisting of at least two species of charged particles: one 
species of ions and one of electrons. 

We also introduce some convenient dimensionless pa- 
rameters which will be used in constructing asymptotic 
orderings for the relevant quantities of the theory. The 
first one, which enters into the construction of the gy- 
rokinetic theory, is defined as Em = max {■ ! -jf-, s — i,e\, 

where r^g — vj_ths/^cs is the species average Larmor 

1/2 

radius, with v± t hs = {T± s /M s } denoting the species 
thermal velocity perpendicular to the magnetic field di- 
rection and tt cs — Z s eB/M s c denoting the species Lar- 
mor frequency. Here L is the characteristic length- 
scale of the inhomogeneities of the EM field, defined 
as L Lb ~ Le, where Lb and Le are the char- 
acteristic lengths of the gradients of the absolute val- 
ues of the magnetic field B (r, t) and the electric field 
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E(r, t), defined as 
-J- = maxll^-lnE 

tures and magnetic fields in AD plasmas, < Em -C 1. 
The second parameter is the inverse aspect ratio de- 



,i = 1,3^. For typical tempera- 



Nested magnetic 
surfaces 
= const 



Central object 



Disc profile 



FIG. 1: Schematic view of the configuration geometry. 



fined as S = ^r^, where Rq is the radial distance from 
the vertical axis to the center of the nested magnetic 
surfaces and r max is the average cross-sectional poloidal 
radius of the largest closed toroidal magnetic surface; see 
Fig.l for a schematic view of the configuration geometry 
and the meaning of the notation introduced here. Then, 
we impose the requirement < <5 <C 1, which is referred 
to as "small inverse aspect ratio ordering". The main 
motivation for introducing this ordering is that we are 
discussing only local solutions where this asymptotic con- 
dition holds; this property also follows from the results 
presented in [l2l EH , and has aready been used in other 
previous work on the subject [l4l Ha ]. The requirement 
S <C 1 is also needed in order to satisfy the constraint con- 
dition imposed by Ampere's law, as discussed in Sec. VII. 
We stress that the S— ordering here introduced is consis- 
tent with the assumption of nested and closed magnetic 
surfaces that are assumed to be localized in space. 

Finally we introduce a parameter St s which measures 
the magnitude of the species temperature anisotropy and 
is defined as St s = l|3 r|| ±3 , where Tj| s and T± s denote 
the parallel and perpendicular temperatures, as mea- 
sured with respect to the magnetic field direction. 

Note that, in the following, we will use a prime " ' " to 
denote a dynamical variable defined at the guiding-center 
position. 



A. The treatment of EM and gravitational fields 

In the following, we shall assume that the magnetic 
field is of the form 



B = V x A = B se// (r, t) + B ext (r,t), 



(1) 



where B se '^ and R ext denote the self-generated magnetic 
field produced by the AD plasma and a non- vanishing ex- 
ternal magnetic field produced by the central object. We 
also impose the following relative ordering between the 
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two components of the total magnetic field: r^psm ~ 

0{e\ lI ), with k > 1. This means that the self-field is 
the dominant component: the magnetic field is primar- 
ily self-generated. Also, the overall magnetic field is as- 
sumed to be slowly varying in time, i.e., to be of the form 
B(r,£ M £), while B sel f and B ext are defined as 

W el f = I(r, e M t)Vif + V^ p (r, s M t) x Vp, (2) 
B ext = V^ D (r,e M t) x Vp, (3) 

where By = I(r, ejat) V<£ and Bp = Vijj p {Y,£Mt) x Vy> 
are the toroidal and poloidal components of the self-field, 
with (tl) p ,(p,'d) defining locally a set of magnetic coordi- 
nates. Moreover, the external magnetic field B ext is as- 
sumed to be purely poloidal and defined in terms of the 
vacuum potential ipD(r,£Mt). In particular, we notice 
here that for typical astrophysical applications of inter- 
est, the function ipu(r, SMt) can be conveniently identi- 
fied with the flux function of a dipolar magnetic field. It 
follows that the magnetic field can also be written in the 
form 

B = /(r,e M t)V(^ + V^(r,e M t) x V<^, (4) 

where the function ip(r,eMt) is defined as i/;(r,ej/t) = 
ipp(v, e m^) + ipD(r,£Mt), and (i/),<p, "&) define a set of 
local magnetic coordinates, as implied by the equation 
B • Vip = which is identically satisfied. In addi- 
tion, it is assumed that the charged particles of the 
plasma are subject to the action of effective EM poten- 
tials {$^(r,£Aft), A(r, £m*)} j where A(r,e M i) is the 
vector potential corresponding to the magnetic field of 
Eq.(gJ), while <& e Jf (r,e M £) is given by 

<P e / f (r,e M t) = $(r,£ M t) + ^^(r.M. ( 5 ) 

Z s e 

with $^(r,£Mt), $(r,£Af*) and 3>g(i*, £m*) denoting 
the effective electrostatic potential and the electrostatic 
and generalized gravitational potentials (the latter, in 
principle, being produced both by the central object and 
the accretion disc). Finally, both the equilibrium effec- 
tive electric field EJ* * , generated by the combined action 
of the effective EM potentials and defined as 

Wjf = -V$f f , (6) 

and the magnetic field B are also assumed to be axi- 
symmetric. 

III. FIRST INTEGRALS OF MOTION AND 
GUIDING-CENTER ADIABATIC INVARIANTS 

In the present formulation, assuming axi-symmetry 
and stationary EM and gravitational fields, the exact 
first integrals of motion can be immediately recovered 
from the symmetry properties of the single charged par- 
ticle Lagrangian function Jzf . In particular, these are the 



total particle energy 

E s = ^ +Zs e^/f(r), (7) 

and the canonical momentum p vs (conjugate to the ig- 
norablc toroidal angle tp) 

Z s e Z s e 
p vs = M s Rv ■ e v H — —ip = — — ip*g. (8) 

Gyrokinetic theory allows one to derive the adiabatic in- 
variants of the system [H, [29[; by construction, these 
are quantities conserved only in an asymptotic sense, 
i.e., only to a prescribed order of accuracy. As is well 
known, gyrokinetic theory is a basic prerequisite for the 
investigation both of kinetic instabilities (see for example 
[32Tl34| ) and of equilibrium flows occurring in magnetized 
plasmas [24. 133-38]. For astrophysical plasmas close to 
compact objects, this generally involves the treatment 
of strong gravitational fields which needs to be based 
on a covariant formulation (see [391-42] ) . However, for 
non-relativistic plasmas (in the sense already discussed) , 
the appropriate formulation can also be directly recov- 
ered via a suitable reformulation of the standard (non- 
relativistic) theory for magnetically confined laboratory 
plasmas [29l - l3ll . I4a449l |. In connection with this, con- 
sider again the Lagrangian function Jzf of charged parti- 
cle dynamics. By performing a gyrokinetic transforma- 
tion of .jSf, accurate to the prescribed order in sm, it fol- 
lows that - by construction - the transformed Lagrangian 
JSf' becomes independent of the guiding-center gyrophase 
angle <j>' '. Therefore, by construction, the canonical mo- 
mentum p'^, s = dJ£"/d<p', as well as the related mag- 
netic moment defined as m' s = jf^P^'si are adiabatic 
invariants. As shown by Kruskal (1962 |5(|) it is al- 
ways possible to determine _£f' so that m' 8 is an adi- 
abatic invariant of arbitrary order in em, in the sense 
that |lnm' s = O + O^/ 1 ), where fl' cs = Z s eB'/M s c 
denotes the Larmor frequency evaluated at the guiding- 
center and the integer n depends on the approximation 
used in the perturbation theory to evaluate m' s . In ad- 
dition, the guiding-center invariants corresponding to E s 
and i/j^s (denoted as E' s and ^ s respectively) can also 
be given in terms of Jz?'. These are also, by definition, 
manifestly independent of <j)' '., 

This basic property of the magnetic moment m' s is es- 
sential in the subsequent developments. Indeed, we shall 
prove that it allows the effects of temperature anisotropy 
to be included in the asymptotic stationary solution. 

Let us now define the concept of gyrokinetic and equi- 
librium KDFs. 

Def. - Gyrokinetic KDF (GK KDF) 

A generic KDF f s (r, v, t) will be referred to as gy- 
rokinetic if its Lagrangian time-derivative 4?f s (r, V, t) 
is independent of the gyrophase angle <f)' evaluated 
the guiding- center position when its state x = (r, v) is 
espressed as a function of an arbitrary gyrokinetic state 
z'= (y', (/)') . More generally, in the following / s (r, v, t) 
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will be referred to as an asymptotic- GK KDF if, neglect- 
ing corrections of order O (e^/ 1 ) ; ^/s( r > v >*) * s inde- 
pendent of (/)' . 

Def. - Equilibrium KDF 

A generic KDF f s (r, v, i) will be referred to as an equi- 
librium KDF if it identically satisfies the Vlavov equation 
-jlfs (r,v,t) = and if f s is also independent of time, 
namely f s = f s (r, v) . More generally, f s (r, v, t) will be 
referred to as an asymptotic- equilibrium KDF if, neglect- 
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fj s (r,v,t) = and 



ing corrections of order O (e 
to the same order f s is independent of t. 

Let us first provide an example of a GK equilibrium 
KDF. This can be obtained by assuming that f s depends 
only on the exact invariants, namely that it is of the 
form f s = fits{E s ,ip*s)i with f* s suitably prescribed 
and strictly positive. On the other hand, an asymp- 
totic GK equilibrium KDF is manifestly of the form 

fs = f*s {E s ,iIj*s,tti's) [again to be assumed as strictly 
positive]. In fact, in this case, by construction, the KDF 
is an adiabatic invariant of prescribed order n, such that 



1 d r 

nL.dt lnf * 



+ (e 



M 



(9) 



(asymptoptic Vlasov equation). In particular, the order 
n (with n > 0) can in principle be selected at will. We 
stress, however, that since gyrokinetic theory is intrin- 
sically asymptotic any GK equilibrium KDF depending 
on the magnetic moment m' s is necessarily asymptotic in 
the sense of the previous two definitions. 

Regarding the notations used in the following, we re- 
mark that, unless differently specified: 1) the symbol "A" 
denotes physical quantities which refer to the treatment 
of anisotropic temperatures; 2) the symbol "*" is used to 
denote variables which depend on the canonical momen- 
tum -i/)* s . 



IV. ASYMPTOTIC EQUILIBRIA WITH 
NON-ISOTROPIC TEMPERATURE 

In this section we derive an equilibrium solution for the 
KDF describing AD plasmas for which the temperature 
is anisotropic. 

Let us first assume that the AD plasma is character- 
ized by a mainly toroidal flow velocity, where the toroidal 
component is expressed in terms of the angular frequency 
by V s (R, z) ■ e v = R,fl s (R, z). In the following, we will 
also require that the plasma is locally characterized by 
a family of nested magnetic surfaces which close inside 
the plasma in such a way that any -0 = const, surface 
is a closed one. We also assume that the system can be 
satisfactorily described by a closed set of fluid equations 
in terms of four moments of the KDF, giving the number 
density, the flow velocity and the parallel and perpendic- 
ular temperatures. 

The presence of a temperature anisotropy means that 
the plasma KDF cannot be a Maxwellian. As already 



mentioned, it remains in principle completely unspeci- 
fied, with just the constraint that it must be a function 
only of the first integrals of motion or the adiabatic in- 
variants of the system. Any non-negative KDF depend- 
ing on the constants of motion and the adiabatic invari- 
ants is therefore an acceptable solution. This freedom in 
choosing a stationary solution is a well-known property 
of the Vlasov equation. We will show that, in these cir- 
cumstances, it is still possible to construct a satisfactory 
asymptotic GK equilibrium KDF (in the sense defined 
in Sec. Ill) which is an adiabatic invariant expressed in 
terms of the two first integrals of motion (JTJ) , © and the 
guiding-center magnetic moment m' s . In particular, the 
form of the stationary KDF which we are going to intro- 
duce is characterized by the following properties: 1) it is 
analytically tractable; 2) it affords an explicit determi- 
nation of the relevant kinetic constraints to be imposed 
on the fluid fields (see the discussion after Eq. (fT5|) ); 3) 
it represents a possible kinetic model which is consistent 
with fluid descriptions of collisionless plasmas character- 
ized by temperature anisotropy; 4) it is suitable for com- 
parisons with previous literature, in which astrophysical 
plasmas have been treated by means of a Maxwellian or a 
bi-Maxwellian KDF (see for example [H, QJ, SI] ) . Then, 
following [l| [24{ , a convenient solution is given by 



(2tt/M, 
x exp 



)3/2 

Hs,, 



T \\*s) 



1/2 



(10) 



(Generalized bi-Maxwellian KDF), where 

r = -22- 

_ &_ 

— — ) 



H = E 



C 



(11) 

(12) 
(13) 



while E s is given by Eq.([7]), ip* 3 is given by Eq.© and 
= t^—- In order for the solution (fTQl) to be 

At s 2j_ s 1 \\* s 

a function of the integrals of motion and the adiabatic 
invariants, the functions /3* s , ct^, Tju s and must 
depend on the constants of motion by themselves. In 
general this would require a functional dependence on 
both the total particle energy and the canonical momen- 
tum. However, in the following, for simplicity, we shall 
consider the case in which only a dependence on ip* s is 
retained [HlHIj]- Namely, /* s depends, by assumption, 



on the flux functions 



rp _ 



a* 



T\\*s (lp*s) ) 

a^s (^*s) , 



(14) 
(15) 
(16) 
(17) 
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which in the following will be referred to as kinetic con- 
straints^ From these considerations it is clear that the 
KDF /* s is itself an adiabatic invariant, and is therefore 
an asymptotic solution of the stationary Vlasov equa- 
tion, whose order of accuracy is uniquely determined by 
the magnetic moment, as already anticipated. 

From definition (|13[) . it follows immediately that an 
equivalent representation for is given by: 



f*s 



P*s exp 



A 



(2ir/M s f 2 (T||„) 



1/2 



J m s (v - v„y 

x exp < — m'.a, 



271 



||*S 



(18) 



where V» s = e ip RQ* s (ip*s) and 



X„ = ( M s + ^Q* 5 Z a e* e Jf ) . (19) 



clear that this asymptotic expansion is valid for r within 
an interval r m ; n < r < r max , where the lower bound is 
fixed by the condition of having a strongly magnetized 
plasma, while the upper bound is given by the geomet- 
ric properties of the system and the small inverse aspect 
ratio ordering. For the purpose of this paper, in per- 
forming the asymptotic expansion we retain the leading- 
order expression for the guiding-center magnetic moment 
,/ _ ,,/ _ m sW ' 2 |gQj]_ Thenj it ig straightforward to 



2B' 



prove that for strongly magnetized plasmas, the following 
relation holds to first order in e (i.e., retaining only lin- 
ear terms in the expansion): /* s = f s [1 + hrj s ] + O (e™), 
n > 2. Here, the zero order distribution f s is expressed 
as 



fs 



(2n/M s f 2 (T lls ) 1/2 T ± 



(20) 



x exp ■ 



M s (v - V s 



2T, 



Ik 



M s w 
2Ai 
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The same kinetic constraints ((Ti)) - pT|) also apply to 
the solution ([T8| . Note that the functions /3* s exp 



M.Rv„ 



-4- 



V* s and T\u a cannot be regarded as fluid fields, since 
they have a dependence on the particle velocity via the 
canonical momentum i))* a . On the other hand, fluid fields 
must be computed as integral moments of the distribu- 
tion function over the particle velocity v. 

Next we show that a convenient asymptotic expan- 
sion for the adiabatic invariant /* s can be properly ob- 
tained in the following suitable limit. Consider, in fact, 
the quantity e defined as e = max {e s , s = i, e}, with 

where we have used the def- 

v> ■ 

inition (JSJ with v v = v • e^, and where L vs denotes the 
species particle angular momentum. We can give an av- 
erage upper limit estimate for the magnitude of £ s in 
terms of the species thermal velocity and the inverse as- 
pect ratio previously defined. To do this, we first set 
ip ~ Bpr 2 , which is appropriate for the domain of closed 
nested magnetic surfaces. Recall that here r is the av- 
erage poloidal radius of a generic nested magnetic sur- 
face. In this evaluation, the species thermal velocities 
Vths and the toroidal flow velocities R£l s are considered to 
be of the same order with respect to the e-expansion, i.e. 
vths/R^s ~ O (e°) (referred to as sonic flow). There- 
fore, assuming v vs ~ v t hs it follows immediately that 
e s ~ where r^ s is the species Larmor radius and 
Lc = rS, with 5 the inverse aspect ratio. We shall 
say that the AD plasma is strongly magnetized when- 
ever < e « 1. This condition is realized if r > r m j n , 

where r m i n = max , s = i, e\ is the minimum aver- 
age poloidal radius of the toroidal nested magnetic sur- 
faces for which e -C 1 is satisfied. In this case e can be 
taken as a small parameter for making a Taylor expan- 
sion of the KDF and its related quantities, by setting 
ip^s ~ ip + O (e fc ), k > 1. From the above discussion, it is 



which we will here call the bi-Maxwellian KDF, where 



and 



the number density n s = rj s exp 



X, 



X„ = M, 



Z £ 

+ — - Z s e<f> 



eff 



(21) 



with rj s denoting the pseudo-density. Then, V s = e v Rrt s 
and the following kinetic constraints are implied from 
d-d}: & = p.ty) = Ty s = T|| S (V), oTs = 



a s (ip) = O s = Qs(ip)- As can be seen, the func- 

tional form of the leading order number density, the 
flow velocity and the temperatures carried by the bi- 
Maxwellian KDF is naturally determined. In partic- 
ular, note that the flow velocity is species-dependent, 
while the related angular frequency fl s must necessar- 
ily be constant on each nested toroidal magnetic surface 
{ip(r) = const.}. Finally, the quantity hrj s represents the 
diamagnetic part of the KDF /* s , given by 



f cM s R 



I Z s e 



Ms_R. 



Y 2 (v-e„) 



(22) 



with Y l = \A U + A 2s - ij - /4 A 4s j , H s = E s - 

^r-ip s £l s (ip s ) and Y 2 = £l s (ip) [1 + ^Ms], where we have 
introduced the following definitions: A\ s = d ^g£° , A 2s = 
A 3s = a^ s = <|L. We remark here 

that: 1) in the e-expansion of (|10[) . performed around 



(|20|) . no magnetic or electric field scale lengths enter, 
as can be seen from Eq. ([22]) : 2) we also implicitly as- 
sume the validity of the ordering -c 1, which will 
be discussed below (see next section). For this reason, 
corrections of Oye^A, with k > 1, to (|2"2"|) have been 
neglected; 3) in this e— expansion we have also assumed 
that the scale-length L is of the same order (with respect 



8 



to e) as the characteristic scale-lengths associated with 
the species pseudo-densities r] s , the temperatures Tju and 
T± s , and the toroidal rotational frequencies il s . 

To conclude this section we point out that the very 
existence of the present asymptotic kinetic equilibrium 
solution and the realizability of the kinetic constraints 
implied by it, must be checked for consistency also with 
the constraints imposed by the Maxwell equations, as 
discussed in Sec. VII. 



V. MOMENTS OF THE KDF 

It is well known that, given a distribution function, 
it is always possible to compute the fluid moments as- 
sociated with it, which are defined through integrals of 
the distribution over the velocity space. Although an 
exact calculation of the fluid moments could be carried 
out (e.g., numerically) for prescribed kinetic closures, in 
this section we want to take advantage of the asymptotic 
expansion of the KDF in the limit of strongly magne- 
tized plasmas to evaluate them analytically, thanks to 
the properties of the bi-Maxwellian KDF. In the follow- 
ing, we provide approximate expressions for the number 
density and the flow velocity, which allow one to write 
the Poisson and Ampere equations for the EM fields in 
a closed form, and for the non-isotropic species pressure 
tensor. Since these fluid fields are then known (in terms 
of suitable kinetic flux functions and with a prescribed 
accuracy), the closure problem characteristic of the fluid 
theories is then naturally solved as well. 

The main feature of this calculation is that the number 
density and flow velocity are computed by performing a 
transformation of all of the guiding-center quantities ap- 
pearing in the asymptotic equilibrium KDF to the actual 
particle position, to leading order in e m (according to the 
order of accuracy of the adiabatic invariant used), and 
they are then determined up to first order in e, in agree- 
ment with the order of expansion previously set for the 
KDF. Terms of higher order, i.e. 0(e^ f ), with n > 1, 
and O (e fe ) , with k > 2, as well as mixed terms of or- 
der O ([££m]") with n > 1, are therefore neglected in the 
present calculation. This approximation clearly holds if 
^r 1 -C 1, which is consistent with the present assump- 
tions. In fact, from the definitions given for these two 
small dimensionless parameters it follows that 

— ~0(<5)«1. (23) 

£ 

To first order in e, the total number density n*°* is 
given by n' ot = J dvf* s ~ n s [l + A„J. Note here 
that the resulting number density has two distinct con- 
tributions: n s is the zero order term given in the pre- 
vious section, while A„ s represents the term of O (e) 
which carries all of the corrections due to the asymp- 
totic expansion of the KDF for strongly magnetized plas- 
mas. The full expression for A„ a is given in Appendix 



A. Finally, a similar integral can be performed to com- 
pute the total flow velocity V* ot . This has the form 
n tot v tot s j dvvjQ n s [V s + AU S ] , where by def- 
inition V s = Q s (Tp)Re v and AU S represents the self- 
consistent FLR velocity corrections given by: 

Aq 

AU S = A^e,, + -^V^ x Vy>, (24) 

where A vs = A„ s Sl s i? + A2 S + A^ s j^. Note that in 
Eq. (l2"4")l the terms proportional to A„ s , A 2s and A 3s 
come from the asymptotic expansion of the KDF for 
strongly magnetized plasmas and are of O (s) with re- 
spect to the toroidal velocity Fl s R. The full expressions 
for A ns , A2 S and A3 S are given in Appendix A. The 
first-order term AU S provides corrections to the zero- 
order toroidal flow velocity with components in all of the 
three space directions and so we can conclude that, al- 
though the dominant fluid velocity is mainly toroidal, 
there is also a poloidal component of order e, associated 
with the term =^-Vip X Vtp. However, this is not nec- 
essarily an accretion velocity, especially under the hy- 
pothesis of closed nested magnetic surfaces which de- 
fine a local domain in which the disc plasma is confined. 
Moreover, note that the ratio between the toroidal and 
poloidal velocities depends also on Sts, hi the sense that 

|v | ~ 0(e) O {Sts)- The magnitude of the tem- 
perature anisotropy can therefore be relevant in further 
decreasing the poloidal velocity in comparison with the 
toroidal one, which on the contrary is not affected by 5t s - 
However, the real importance of this result in connection 
with the astrophysics of collisionless AD plasmas is, in- 
stead, the fact that this poloidal velocity is a primary 
source for a poloidal current density which in turn can 
generate a finite toroidal magnetic field (see the section 
on the Maxwell equations). This means that, even with- 
out any net accretion of disc material (which would re- 
quire at least a redistribution of the angular momentum) , 
the kinetic equilibrium solution provides a mechanism for 
the generation of a toroidal magnetic field, with serious 
implications for the stability analysis of these equilib- 
ria. The physical mechanism responsible for this poloidal 
drift is purely kinetic and is essentially due to the conser- 
vation of the canonical toroidal momentum and the FLR 
effects associated with the temperature anisotropy. As 
a last point, consider the species non-isotropic pressure 
tensor, which is defined by the following moment of the 
KDF: g s = fdvM a (v - V*°*)(v - V*°*)/T s . Then, the 
overall pressure tensor of the system is obtained by sum- 
ming the single species pressure tensors: II = J2 s =i eH s - 
A direct analytical calculation retaining only the zero- 
order terms (with respect to all of the small dimension- 
less parameters) in the Taylor expansion of the KDF /* s 
(whose leading-order expression coincides with the bi- 
Maxwellian KDF), shows that the corresponding species 
tensor pressure is non-isotropic to leading order and in 
this approximation is given by: 

E s =p±sl+(P\\ S -P± s )bb, (25) 



9 



where p± s = n s Tj_ s and p|| s = n s T\\ s represent the 
leading-order perpendicular and parallel pressures. The 
divergence of the species pressure tensor is of particular 
interest; this is given by: 

V • g s = V P±S + bB • V ( P " a ~ P±3 ) - Ap a Q, (26) 

where Q = [bb- V In B + i§b x J - V In B] and Ap s = 

(P\\s ~P±s)- 



VI. THE CASE OF ISOTROPIC 
TEMPERATURE 

In this section, we consider the case of isotropic tem- 
perature for the equilibrium distribution /* s . When the 
condition Tm s = T± s = T s is satisfied, the stationary 
KDF reduces to /* s , where 

/**= .^ expl-^j (27) 



(2T* S /M S ) 3/2 



T 

: ! : 



is referred to as the Generalized Maxwellian Distribution 
with isotropic temperature [l5| . Here, H* s retains its defi- 
nition (|13[) . while the kinetic constraints are expressed for 
the quantities n* s and T* s , whose functional dependence 
is ?7*s = ri s {ip* s ) and T* s = T s (ip* s ). By construction, 
this distribution function is expressed only in terms of 
the first integrals of motion of the system and is there- 
fore an exact kinetic equilibrium solution. Performing 
an asymptotic expansion in the limit of strong magnetic 
field, as done before for /* s , gives the following result: 
Us = fit. [1 + h Ds ] + 0(e n ), n>2, where 

( M s (v-V s ) n 



fhls — 



7T 3 / 2 (2T S /M 8 ) 3/2 



exp • 



2T, 



(28) 



is the zero-order term of the series, which coincides with 



a drifted Maxwellian KDF with T, 



Q s (?p)Re v and n s = rj s (ip) exp 
function ho s is given by 



In this case, the 



cM s R 
Z s e 



M S R 



Y, 



with Y\ = 

si.(i/>) [l + Mss] 

^ _ ainn.(^) 



where A\ t 



dip 



Ms. _ 3 
T s 2 

d In rj 

= — £hp~ 

Finally, as shown in [l 



and 

A 2s 



(29) 
Y 2 = 

_ In T 3 
— dip ' 

the angu- 

9<x: 



lar frequency is given to leading order by fl s (tp) — -g^j 

where % = c&'jff + |£| lnn s . 

Before concluding this section, we stress again that the 
case with isotropic temperature represents a result whose 
accuracy is not limited by dependence on any gyrokinetic 
invariant and that it does not require any guiding-center 
variable transformation. For this reason, there are no 
restrictions of applicability of the solution (|27|) which, in 
principle, holds also in the limit of vanishing magnetic 
field. 



VII. THE MAXWELL EQUATIONS AND THE 
"KINETIC DYNAMO" 

In this section we write the Poisson and Ampere equa- 
tions for the EM fields explicitly, pointing out the con- 
sequences of the kinetic treatment developed in Sec- 
tions IV- VI. In particular we prove that, besides a self- 
generated poloidal magnetic field, the kinetic equilibrium 
can sustain also a toroidal field (which may be thought 
of as being a kinetic dynamo), thanks to the combined 
effects of FLR corrections and temperature anisotropics. 
For definiteness, let us consider the Poisson equation for 
the electrostatic potential <f>, expressed as 



V 2 $ 



-4-7T 



E 



q s n t 



[1 + A r 



(30) 



where A„ s is written out explicitly in Appendix A. In 
the limit of a strongly magnetized plasma and consid- 
ering the accuracy of the previous asymptotic analyti- 
cal expansions, we shall say that the plasma is quasi- 



neutral if the ordering 



+ O(e k 



with k > 2 holds, whereas we call it weakly non-neutral 



if 



= + O (e) . The kinetic equilib- 



rium for a weakly non-neutral plasma is referred to as a 
Hall kinetic equilibrium [l5[ , and the corresponding fluid 
configuration is referred to as a Hall Gravitational MHD 
(Hall-GMHD) fluid equilibrium [13 ]. 

Next, we show that quasi-neutrality (in the sense just 
defined) can be locally satisfied by imposing a suitable 
constraint on the electrostatic (ES) potential $. It can 
be shown that this constraint can always be satisfied 
since the leading-order contribution to the ES poten- 
tial remains unaffected. The result follows by neglecting 
higher-order corrections to the number density (A„J and 
setting qi — Ze and q e — —e. Thanks to the arbitrari- 
ness in the choice of the flux functions introduced by the 
kinetic constraints (see Sec. IV), it is possible to show 
that quasi-neutrality implies the following constraint for 
the oscillatory part $~ of the ES potential, i.e., correct 
to bothO(e°) andO(e^), 



(31) 



where S" 



Xj 



and X s 



the constraint (S~) 
follows that 

Vs 



ipQ, s (tfi) — M s $>g ) . In particular, the ar- 

Zr,i 



bitrarincss in the coefficient Szr can be used to satisfy 



0. In fact, in view of Eq. fT^]) . it 

Ps MT|| a (V0 



1 



(32) 



where a s (ip) is related to [ijj) as outlined in Appendix 
A and the flux functions still remain arbitrary. In conclu- 
sion, Eq. (|3"Tl) determines only $~ and not the total ES 
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potential. Note that this solution for the electrostatic 
potential $~ can be shown to be consistent with earlier 
treatments appropriate for Tokamak plasma equilibria 
[24l l36j . This can be exactly recovered thanks to the ar- 
bitrariness in defining the pseudo-densities and by taking 
the limit of isotropic temperatures and zero gravitational 
potential, as in the case of laboratory plasmas. In this 
limit the species pseudo-densities become flux functions 
[24| . Then, because of this arbitrariness, by taking 



Ve 

Zrji 



1, 



(33) 



it follows that Eq. (f3"Tj) reduces to the form presented in 
)24j, which can only be used to determine the poloidal 
variation of the potential. 

Let us now consider the Ampere equation. Adopting 
the Taylor analytical expansion of the asymptotic equi- 
librium KDF and neglecting corrections of O (e k ) , with 
k > 2, and O (e^ 1 ), with n > 0, this can be approxi- 
mately written as follows: 

4-7T 

s—i,e 



- AU S ] , (34) 



where B se '^ is as defined in Eq.Q and the expression for 
AU S is given by Eq. (j2~4"|) . The toroidal component of this 
equation gives the generalized Grad-Shafranov equation 
for the poloidal flux function tp p : 

A*^ p = -—R^q s n s [n s {iP)R + A vs ], (35) 



where the elliptic operator A* is defined as A* = R 2 V ■ 
(i? _2 V) [5l|. The remaining terms in Eg. (I3~4l give the 
equation for the toroidal component of the magnetic field 

7(^,i?) 



R 



. In the same approximation, this is: 



V/^, J?) x 



4tt 



,n.-5 i VV» x Vtp, (36) 

73 



field without requiring any net accretion and in the ab- 
sence of any possible instability /turbulence phenomena. 
This new mechanism results from poloidal currents aris- 
ing due to the FLR effects and temperature anisotropics 
which are characteristic of the equilibrium KDF. We re- 
mark that the self-generation of the stationary magnetic 
field is purely diamagnetic. In particular, the toroidal 
component is associated with the drifts of the plasma 
away from the flux surfaces. In the present formula- 
tion, possible dissipative phenomena leading to a non- 
stationary self field have been ignored. Such dissipative 
phenomena probably do arise in practice and could oc- 
cur both in the local domain where the equilibrium mag- 
netic surfaces are closed and nested and elsewhere. Tem- 
perature anisotropies are therefore an important physi- 
cal property of collisionless AD plasmas, giving a possi- 
ble mechanism for producing a stationary toroidal mag- 
netic field. We stress that this effect disappears alto- 
gether in the case of isotropic temperatures, as demon- 
strated in Appendix A. Finally, we consider the ratio 
between the toroidal and poloidal current densities (Jt 
and Jp). In the small inverse aspect ratio ordering, ne- 
glecting corrections of order O (<5 fc ) with k > 1, this pro- 
vides an estimate of the magnitude of the corresponding 
components of the magnetic field. In fact, in this limit 



we can write 



VxB r | 



|B T | 



|Jf 



and so conclude 



|VxB P ]W^] TJ7T 
that, although for the single species velocity, the order- 

^V-i/jx V<p\ 



mg 



O (e) O {S Ts ) holds (see Sec. VI), this 
might no longer be the case for the magnetic field, which 
instead depends on the ratio between the total toroidal 
and poloidal current densities. In particular, the possi- 
bility of having finite stationary toroidal magnetic fields 
is, in principle, allowed by the present analysis, depend- 
ing on the properties of the overall solution describing 
the system. 



VIII. CONCLUSIONS 



where A 3s , given in Appendix A, contains the contribu- 
tions of the species temperature anisotropies. For con- 
sistency with the approximation introduced, in the small 
inverse aspect ratio ordering, it follows that aJ g^ ,1? ^ = 
+ O (6 k ) , i.e., to leading order in 5: I = I(ip) + (6 k ) , 
with k > 1. This in turn also requires that the corre- 
sponding current density in Eq.(|36|) is necessarily a flux 
function. Then, correct to 0(e), 0(eP m ) and O (<5°) , the 
differential equation for / (ip) becomes: 



dl(i>) _ 4tt 
dip c 



q s n s — 



(37) 



which uniquely determines an approximate solution for 
the toroidal magnetic field. This result is remarkable be- 
cause it shows that there is a stationary "kinetic dynamo 
effect" which generates an equilibrium toroidal magnetic 



Getting a complete understanding of the dynamical 
properties of astrophysical accretion discs still represents 
a challenging task and there are many open problems 
remaining to be solved before one can get a full and con- 
sistent theoretical formulation for the physical processes 
involved. 

The present investigation provides some important 
new results for understanding the equilibrium proper- 
ties of accretion discs, obtained within the framework 
of a kinetic approach based on the Vlasov-Maxwell de- 
scription. The derivation presented applies for collision- 
less non-relativistic and axi-symmetric AD plasmas un- 
der the influence of both gravitational and EM fields. 
A wide range of astrophysical scenarios can be investi- 
gated with the present theory, thanks to the possibility 
of properly setting the different parameters which char- 
acterize the physical and geometrical properties of the 
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model. A possible astrophysical context is provided, for 
example, by radiatively inefficient accretion flows onto 
black holes, where the accreting material is thought to 
consist of a plasma of collisionless ions and electrons with 
different temperatures, in which the dominant magnetic 
field is generated by the plasma current density. We have 
considered here the specific case in which the structure 
of the magnetic field is locally characterized by a fam- 
ily of closed nested magnetic surfaces within which the 
plasma has mainly toroidal flow velocity For this, we 
have proved that a kinetic equilibrium exists and can 
be described by a stationary KDF expressed in terms 
of the exact integrals of motion and the magnetic mo- 
ment prescribed by the gyrokinetic theory, which is an 
adiabatic invariant. Many interesting new results have 
been pointed out; the most relevant ones for astrophys- 
ical applications are the following: 1) the possibility of 
including the effects of a non-isotropic temperature in the 
stationary KDF; 2) the proof that the Maxwellian and bi- 
Maxwellian KDFs are asymptotic stationary solutions, 
i.e. they can be regarded as approximate equilibrium 
solutions in the limit of strongly magnetized plasmas; 3) 
the possibility of computing the stationary fluid moments 
to the desired order of accuracy in terms of suitably pre- 
scribed flux functions; 4) the proof that a toroidal mag- 
netic field can be generated in a stationary configuration 
even in the absence of any net accretion flow if and only 
if the plasma has a temperature anisotropy. 

This last point, in particular, is of great interest be- 
cause it gives a mechanism for generating a stationary 
toroidal field in the disc, independent of instabilities re- 
lated to the accretion flow. The consistent kinetic for- 
mulation developed here permits the self-generation of 
such a field by the plasma itself, associated with localized 
poloidal drift-currents on the nested magnetic surfaces as 
a consequence of temperature anisotropics. This station- 
ary poloidal motion is made possible in the framework of 
kinetic theory by the conservation of the canonical mo- 
mentum as a result of FLR effects. 

These results reveal and confirm the power of the ki- 
netic treatment and the necessity for adopting such a 
formalism in order to correctly understand the physical 
phenomena occurring in accretion discs. This study may 
represent a significant step forward for understanding 
the physical properties of accretion discs in their kinetic 
equilibrium configurations, and it motivates making fur- 
ther investigations of the subject aimed at extending the 
present range of validity to more general physical config- 
urations. The conclusions presented here may also have 
important consequences for other applications, and can 
provide a reference starting point for future work on ki- 
netic stability analysis of accretion discs. 
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Appendix A: Calculation of the fluid moments 

In this appendix we give the detailed expressions for 
the coefficients which determine the fluid moments of the 
KDF, obtained by integrating the KDF over the velocity 
space. For the number density, we have found that n s ot = 
n s [1 + A„J. The term A„ s is given by 



A„ s 



where V s 



V, 



7i + 73 



Ti 



4T 



±8 



2 73 / 2 (T|, s - T 



M s 

± .i) V s 



v: 



(Al) 



7i 



72 



73 



B 2 R 2 M S 

M7 S (V>) and 

cM s R M s \ 
-K 



Z s e 
cM s R 

Z s e 
cM 2 R A 2 



\\s 



A 4i 



Z,e 2Ti 



in which K 

da s 



A ls + A 2s 



■ft. 



(A2) 
(A3) 
(A4) 

and 



Ais = , with a s (ip) = 2f~-. Note that here a s (i/i) 
differs from a~g (ip) because of the guiding-center trans- 
formation of the magnetic field B. 

Proceeding in the same way for evaluating the second 
moment of the KDF, it can be shown that the first order 
correction AU S to the toroidal flow velocity can be writ- 



ten as AU S = A^e^ 
that 



x V(y5, where we recall 



A vs = A ns n s R + A 2s 



A ^RB- ^ 

Here A„ s is as given in Eq. (|Aip , while A2 S and A3 S are 
given by 



A 2s = 



^(7i + 3 7 3K 2 )- 



BAL 



A 



3 s 



RB 2 M t 



(2T ±S 

- t, ; 



4T ±S ) 
— T\i.) 



(A6) 



(A7) 



RBM. 
I73 



(71 + 373 K 2 



RBM 2 V 3T " S 



ATj 



■ T\\ S T± S 
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Finally, in the limit of isotropic temperatures, the solu- 
tion for the number density is 



where now 72 = 0, 
7i = 



f cMgR A 2 
73 = 1 Z s e 2T a f' 



and K = 



A-ls + J\2s \ ji 2 



1) 



(A8) 



cM s R r _ M S V S .„ , , . . 
— — ^ + [1 + , (A9) 

Z s e l s 



(A10) 



. For 



calculating the flow velocity from Eq. (|A5l) . in the same 
limit, A„ s is as given by Eq. (|A8p . A 2s reduces to 



(All) 



and A 3s = since, in Eq. (|A7p . the second and third 
terms on the right hand side necessarily vanish, while the 
first one is proportional to = 7^7% where a s (ip) = 



b _ 



A, 



= 0, and hence vanishes too. 
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